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Abstract 

The superconducting state of underdoped cuprates is often described in terms of a single energy- 
scale, associated with the maximum of the (d-wave) gap. Here, we report on electronic Raman 
scattering results, which show that the gap function in the underdoped regime is characterized 
by two energy scales, depending on doping in opposite manners. Their ratios to the maximum 
critical temperature are found to be universal in cuprates. Our experimental results also reveal 
two different quasiparticle dynamics in the underdoped superconducting state, associated with two 
regions of momentum space: nodal regions near the zeros of the gap and antinodal regions. While 
antinodal quasiparticles quickly loose coherence as doping is reduced, coherent nodal quasiparticles 
persist down to low doping levels. A theoretical analysis using a new sum-rule allows us to relate the 
low-frequency-dependence of the Raman response to the temperature-dependence of the superfluid 
density, both controlled by nodal excitations. 

PACS numbers: 74.72.-h, 74.62.Dh, 78.30.-j 
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Understanding the origin of high temperature superconductivity in copper oxide based 
materials is one of the outstanding problems in modern condensed matter theory. Two 
decades of intensive theoretical and experimental work have revealed that this phenomenon 
takes place in various families of cuprates, all containing copper oxide layers separated 
by building blocks that provide a reservoir of carriers. The cuprates phase diagram is 
remarkably universal. When there is nominally one hole per copper in the copper oxygen 
planes the cuprates are in an insulating phase due to correlations known as a Mott insulator. 
As the number of charge carriers is increased, a process called doping, the cuprates turn 
into a d-wave superconductor. 

The superconducting temperature has a dome like shape as a function of doping and 
cuprates exhibit two physically distinct regimes: i) the under-doped regime where the criti- 
cal temperature Tc increases with doping, ii) an over-doped regime where Tc decreases with 
doping. They are separated by optimal doping where Tc reaches its maximum. There is 
now consensus that the "normal state" above the superconducting critical temperature in 
the under-doped and optimally doped regimes is far from being normal. In particular, a 
pseudo gap develops in the underdoped regime, corresponding to a partial suppression of 
spin and charge excitations. This phenomenon lies outside the standard theory of solids 
and manifests itself in all experiments, such as neutron scattering cross sections, nuclear 
magnetic resonance (NMR), magnetic susceptibility, specific heat, angle resolved photoe- 
mission spectroscopy (ARPES), optical conductivity, tunneling and Raman spectroscopies 
(for a review, see Ref. 

In contrast, there is no consensus on whether the superconducting state in the underdoped 
regime, which emerges from the Mott insulating state, can be described by the standard 
BCS-Migdal-Eliashberg theory for a d-wave superconductor. In this article, we address this 
issue using Electronic Raman Spectroscopy (ERS). ERS is a powerful tool for probing quasi- 
particles of the superconducting state in selected parts of the momentum space namely the 
antinodal (ANR) and nodal (NR) regions where the amplitude of the superconducting gap 
reaches its maximum and vanishes respectively 0. Here ERS is performed on HgBa2Cu04+5 
(Hg-1201) one of the most simple structure which takes a pure tetragonal symmetry with 
only one copper oxide plane per unit cell. By a judicious choice of the laser excitation lines 
we have enhanced the ERS signal allowing a direct view of the electronic properties of this 
cuprate without invoking ad hoc phonon subtraction procedures usually used in cuprate 
systems. 

Our results demonstrate that the superconducting state on the underdoped side is also 
not normal, and cannot be described within the standard d-wave BCS model. We show 
that a minimal description of this state requires two energy scales and two quasiparticle 
dynamics. One scale is associated with the nodal region, and decreases with decreasing dop- 
ing, while the second is associated with the antinodal region, and increases with decreasing 
doping. Close to the nodes the quasiparticle spectral weight remains substantial even at 
low doping in contrast to the antinodes where it drastically decreases. As a consequence 
i) we show that the superconducting order parameter is controlled by two parameters (the 
slope of the gap function at the nodes and the antinodal gap maximum) instead of one 
as expected from a standard d-wave gap, ii) we find a strong momentum dependence of 
the quasiparticles spectral weight in the superconducting state of the underdoped regime, 
and iii) we finally establish, using a new Raman sum-rule, a simple relationship between 
the superfluid density and the low-energy Raman scattering associated with nodal physics, 
suggesting that the Fermi liquid renormalization of the current (vector quantity) and stress 
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tensor (which transforms as a product of currents) have similar doping dependence. Our 
new experimental results place strong constraints on current theories of the high temperature 
superconductivity phenomenon. 

EXPERIMENTAL RESULTS 

The ERS measurements (see Methods section for details) have been performed on se- 
lected, as-grown, Hg-1201 single crystals from several batches with different doping lev- 
els extending from the slightly over-doped (OD) to the underdoped (UD) regime: (Tc = 
92K(OD), 95K(OP), 89K(UD), 78K(UD) and 63K(UD)). 

ERS is an energy probe, but also a momentum probe which allows us to select different 
parts of the Fermi surface by choosing combinations of incident and scattered light polar- 
izations 0. The NR and ANR have been explored using cross polarizations parallel to the 
Cu-0 bond directions (B23) and at 45° from them (Big), respectively. All spectra have been 
corrected for the spectral response of the spectrometer and for the Bose-Einstein factor. 
They are thus proportional to the imaginary part of the Raman response function x"(a;). 

Figure Q displays the NR (B2g) and ANR (Bi^) Raman response functions in both the 
normal and superconducting states of the Hg-1201 single crystals, at various doping levels. 
At optimal doping (Tc=95 K), the electronic Raman continua for both the ANR and NR 
exhibit a redistribution of spectral weight from energies lower than 400 cm~^ to higher 
energy, when going from the normal state to the superconducting state. This redistribution 
is more pronounced for the ANR than for the NR. At low energy (below 400 cm~^), the 
ANR superconducting continuum exhibits a cubic frequency dependence with a well-marked 
superconducting pair breaking peak, indicated by an arrow on FigUl at a frequency Huan — 
505 cm~^ (~ SksTc). In contrast, the superconducting spectrum in the NR displays a linear 
frequency dependence up to 400 cm~^, as well as a weaker signature of the pair-breaking peak 
(arrow) close to the same frequency /icjat ~ Uujan than for the ANR spectrum (but with a half 
width at half maximum reduced approximately by a factor of two). The Ramanresponse 
at optimal doping is thus characterized by a single energy scale fioJAN — Tiojn associated 
with the pair-breaking peak, and all the features described above are consistent with those 
expected for a d-wave superconductor with a maximum value of the superconducting 
gap given by 2Am = ^ujan (see below). Our results for one overdoped sample (spectra at 
the top of Fig. Q]) can also be interpreted in terms of a single energy scale. 

In contrast, as doping is decreased below the optimal one, the evolution of the Raman 
spectra in the superconducting phase becomes strikingly different in the ANR and in the 
NR. As the doping level (and Tc) is reduced, the characteristic energy of the antinodal 
peak (indicated by an arrow on the right pannel of Fig. increases. Simultaneously, the 
intensity of this peak rapidly decreases as Tc decreases, and finally disappears in the vicinity 
of Tc=78K. In contrast, the characteristic energy scale of the spectrum in the NR (which we 
take to be the frequency oi the maximum observed in 'x"{uj)) follows T^.. Furthermore, 
contrarily to the antinodal peak, the nodal peak persists down to the lowest doping that we 
have studied (Tc=63K). We note that similar observations have been reported previously 
on other cuprate materials, such as Bi2Sr2CaCu208+5 (Bi-2212), Bi2Sr2_a;Laa;Cu06+5 (Bi- 
2201), YBa2Cu307-5 (Y-123) and La2-^Sr^Cu04 (LSCO) 0, Q. This demonstrates that 
the electronic Raman response in the underdoped regime involves two distinct energy scales, 
with opposite doping dependence. As discussed below, this is inconsistent with a simple 
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BCS d-wave description [3|. 

In order to substantiate further this point, we have plotted in Fig. |21the characteristic 
ratios huAj^j /ksT"^"'^ and huj]^ /ksT^""^ obtained for several different families of cuprates by 
different groups as a function of doping at a fixed temperature well below T^. The 

doping value p is inferred from Tc using Tallon's equation 0: = 82.6 (p-0. 16)2, 

with T'^"'^ the critical temperature at optimal doping. Fig. |21 reveals that these ratios have a 
universal dependence on doping. For underdoped compounds, two distinct scales are present, 
with the two ratios behaving in opposite manners as a function of doping, while a unique 
energy scale and doping dependence is recovered at optimal doping and in the overdoped 
regime. 



INCONSISTENCY WITH A SIMPLE BCS MODEL 

Let us now analyze these results using the simplest possible framework, that of a BCS 
superconductor with a (i-wave gap function of the form = cos(2</)) (where (f) is the 
angle ^associated with momentum k on the Fermi surface). The Raman response would then 
read 



_ 2vriV^ / (7^^-^(0))^A^cos2(20) \ 



^ \ v/^2_4A2^eos2(20) I 

7'^^'^(0) is the Raman vertex associated with each polarization: 7"^^(0) = 'ysig cos(20) 
while 7^(0) = 7B2g sin(20) , and {■..)pg denotes a Fermi-surface average. This predicts 
a sharp pair-breaking peak (corresponding to a divergence of this expression) in the Big 
(ANR) geometry at a; = 2Am, and a weaker singularity in the B2g (NR) geometry at the 
same frequency scale (furthermore, within the simplest BCS formula above, the i?2g (NR) 
response displays a maximum at a somewhat lower energy than the peak in the Big (ANR) 
channel, but both are governed by one energy scale, that of the maximal gap A^). This 
is roughly consistent with the experimental observations at optimal doping, but completely 
fails to account for the underdoped spectra (for which the observed peaks become clearly 
distinct in each polarization and have opposite doping dependence). 

This clearly demonstrates that one or both of the following assumptions become invalid 
in the underdoped regime: i) non-interacting BCS quasiparticles ii) a gap function with 
the simple form A^ = AmCos(2</)) characterized by a single energy scale. Moving away 
from assumption i) requires taking into account, in the framework of the Landau theory of 
interacting quasiparticles, the spectral weight of these quasiparticles, in general smaller 
than one and k-dependent, as well as the Fermi-liquid vertex A^ describing the interaction 
of the quasiparticles with external perturbations. This leads to: 

in which a general gap-function has also been taken into account. This expression contains 
two unknown functions of momentum on the Fermi-surface however: (ZA)^ and A^, and 
constructing on this basis a phenomenological analysis (to which we shall come back later 
in this paper) requires further insight. 
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LOW-ENERGY NODAL EXCITATIONS, SUM- RULE AND RELATIONSHIP 
TO THE SUPERFLUID DENSITY 



In order to gain such insight, we have focused on the low-energy part of the Raman 
spectra, which is controlled by the properties of the nodal quasiparticles. The 32^, geometry 
is particularly significant in this respect, since the largest low-energy response is obtained 
in this geometry and the NR is directly sampled. Figures Q and IHl demonstrate that a linear 
dependence on frequency is found in this geometry, for several dependent doping levels. This 
is expected from the above formula, which yields: Xn{^ — > 0) = 7^23 i^-^)!^^ + •••• 
this expression, = is the slope of the gap function at the nodes, and {ZA)n is the 
value of {ZA)k at the node. Hence, a study of the doping- dependence of the slope of the nodal 
(B2g) response, allows determining, in principle, the important parameter a = ^(ZA)^ 
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associated with nodal physics. 

In order to compare the slopes of samples with different doping levels however, the Raman 
spectra must be properly normalized. To achieve this goal, we have established a theoretical 
sum-rule for the Raman intensity of a weakly doped Mott insulator, which reads, at low 
doping levels p: x"{^)^duj = Cp. We have derived this expression (see Methods section 
below) by starting from a one-band t-J model, in which case the upper cutoff Q can be taken 
to infinity, and the constant C depends on the hopping amplitudes and exchange constant: 
C{t,t',J). We expect this expression to be general for a weakly doped Mott insulator, 
provided the upper cutoff is taken to be of the order of the bare bandwidth and significantly 
lower than the scale U associated with the upper band. We have normalized the Raman 
data using this sum-rule for our own data on Hg-1201. This is shown in Fig. El and reveals 
that in the doping range of interest the slope is doping- independent. In order to point out 
the universality of this phenomena, we have applied the same Raman sum-rule to a number 
of spectra previously published in the literature for Y-123 and Bi-2212 The results 

for the slope a as a function of doping are displayed on Fig. |31 It is seen that the slope is 
essentially doping- independent for all doping levels between optimal doping (p=0.16) down 
to p=0.09. For smaller doping levels, the Y-123 data (two data points available) suggest 
that a may end up decreasing at very low doping. 

This doping-independence of the Raman (NR) slope over an extended range of doping 
levels is very reminiscent of the behavior of the slope of the linear term in the temperature- 



dependence of the penetration depth (or superfluid density) |10|, |l]|: Ps{T) = ps{0) — 
[3T This quantity is also associated with the physics of nodal quasiparticles, and 

given by almost the same expression [l^, E3] oc ^(ZAp)^, except for the fact that 
the Fermi liquid parameter Ap corresponds to a different angular-momentum channel than 
in the case of Raman. The comparison made in Fig. |3] between these two quantities shows 
excellent agreement, hence establishing a previously unforeseen relation between Raman and 
penetration depth measurements. We also note that very recent experiments 0| do suggest 
deviations from (3 = const, at very low doping levels, in agreement with the two highly 
underdoped Y-123 data points in Fig. HJ As pointed out in Ref. l3|, the independence of 
the slope P of the superfluid density upon doping, in an extended range, is a serious difficulty 
for the simplest ("vanilla") version of RVB theory in either its auxiliary boson or 
wave function formulations ^3] observations reveal a similar problem in connection 

with Raman scattering. 

Having established the doping- (in) dependence of the a = ^(ZA)^ parameter associated 
with nodal physics, we can make further progress in disentangling the relative effects of i) 
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quasiparticle dynamics and ii) the form of the gap function on the Raman spectra. Let us 
first assume that the gap function has the simple form: = cos(20) characterized by 
a single energy scale (the maximum gap), and such that = SA^. Since the behavior of 
the ANR (Big) peak implies that Am increases as doping is reduced (Fig. the doping- 
independence of a would then imply that {ZA)n must correspondingly increase as p is 
reduced. This behavior is highly unlikely to apply for a doped Mott insulator, since it would 
correspond to a reinforcement of the quasiparticle spectral weight as doping is reduced. 
Furthermore, ARPES data and theory do suggest that the nodal quasiparticle spectral 
weight decreases as p is reduced [l^ . Hence, we are led to suspect that the slope of the 
gap at the nodal points, , must decrease as the doping level is reduced (in order to keep a 
constant) and hence does not track the maximum gap A.^. This, of course, has far reaching 
consequences, namely: that a pure cos (20) form does not hold in the underdoped regime, and 
that two energy scales characterize this regime (as already foreseen from the above analvsis 



in Fig. EJ. Indeed, this has been previously suggested from ARPES experiments 21, 
indicating that the superconducting gap function may change from a "V-shape" to a "U- 
shape" as the doping level is reduced. 

PHENOMENOLOGICAL DESCRIPTION 

We now return to the general expression for the Raman response based on Fermi liq- 
uid considerations, and explore whether reasonable momentum- and doping-dependence 
of the two functions ZA{(f)) and A(0) (associated respectively with quasiparticle physics 
and with the superconducting gap) can describe our data. We use the shape of the gap 



function (consistent with (i-wave symmetry) proposed from ARPES data in Ref. j21 
A(0) = Am(i?cos20 + (1 — i?)cos60), where < i? < 1. This form is characterized 
by two scales: the (antinodal) maximum gap A^ and the nodal slope f a = Am(8i? — 6). In 
a first attempt, we have taken ZA to be ^-independent (uniform along the Fermi surface). 
The parameters A^ and B were varied, as in Ref. jllj, and the values of ZA chosen such as 
to keep the ratio {ZA)j^/vA constant, for consistency with our experimental observation on 
the doping- independence of the slope a. We found that the resulting spectra (not shown) 
could capture only part of our experimental observations: the Bi^ peak does move to higher 
energy but does not loose intensity quickly enough, and the response becomes very 
flat as V/^/Am is reduced, but its maximum does not clearly recede towards lower energy. 
Hence, a variation of ZA{(j)) along the Fermi surface is clearly needed, which we have chosen 
as displayed in Fig. (chosen again in such a way that (ZA)^/fA is kept constant). The 
four curves colors in Fig. El are meant to correspond, qualitatively, to the overdoped (BCS), 
optimally doped, slightly and strongly underdoped regimes. The resulting spectra are dis- 
played in Figs.Efc-d for both Big and B2g geometries. All the main qualitative aspects of our 
experimental observations are captured by this simple phenomenological model, namely: the 
opposite variation of the nodal and antinodal peaks with energy, the constant nodal slope, 
and the drastic suppression of the antinodal peak due to the incoherence of the antinodal 
quasiparticles (consistent with Fermi surface "arcs" (23l|). 
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COMPARISON TO OTHER EXPERIMENTS 



Finally, we discuss whether other experimental probes support the existence of two energy 
scales in the superconducting state of the underdoped cuprates. 

Let us first note that the two basic building blocks of the phenomenological analysis 
developed above are consistent with available information from ARPES experiments. The 
leading edge in the ANR (corresponding to the higher energy scale) was shown to increase as 
the doping level is reduced ji^l • Evidence for a second energy scale associated with the nodal 
region, and varying in an opposite manner with doping, has been suggested from ARPES 
in Refs. j^, 2^, with a corresponding change of the gap function as discussed above. More 
extensive ARPES studies of the leading edge in the NR are clearly needed. ARPES also 
reveals that the (gapped) antinodal quasiparticles in the superconducting state quickly loose 
spectral weight and eventually become incoherent as doping is reduced (20, |2^, while nodal 
quasiparticles loose s pec tral weight as doping is reduced but maintain coherence down to 
low doping levels [ifll l2fll|. 

It was suggested from thermal conductivity measurements on Y-123 jl^ that the nodal 
slope of the gap, increases as the doping level is reduced, which would seem in contra- 
diction with our observations. However, closer examination [27[ reveals that only one data 
point is available in the doping range of interest here, and more experiments are clearly 
needed. 

Tunneling, while not being a momentum- selective probe like Raman or ARPES, exhibits 
a set of distinct superconducting line shapes as a function of doping. The main peak in 
these spectra moves to higher energies as doping is reduced from optimal doping |28L lioj ] 
(as also seen from break-j unction measurements jl^ 0]), and at the same time broadens 
considerably. This is consistent with the behavior of the larger energy scale, and with 
the rapid loss of coherence of antinodal quasiparticles. Fourier transform tunneling has 
shown that antinodal decoherence is closely related to a local charge order, at least in 
surface ji^. This is consistent with the loss of the electronic Raman response which is 
sensitive to the charge dynamics. In contrast, the observation of quasiparticle interferences 
from tunneling provides strong evidence for the coherence of nodal quasiparticles. There 
are also some recent indications that the lower energy scale might show up as a low-voltage 
shoulder in the tunneling spectra of underdoped samples [l^ . There are also some previous 
indications in the literature that the lower energy scale shows up as the gap edge in the 
Andreev reflection spectra of underdoped samples js^]. 



CONCLUSION 



In summary, electronic Raman scattering clearly demonstrates the existence of two dis- 
tinct energy scales in the superconducting state of underdoped cuprates, depending on dop- 
ing in opposite ways. Here, we have suggested that these two scales are associated with the 
nodal slope of the gap function, and its maximum, respectively. Correspondingly, our exper- 
iments reveal two different dynamical properties of the quasiparticles in these two regions. 
Hence, the dichotomy between the coherence of nodal quasiparticles and the incoherence of 
antinodal ones, which is well established in the normal state, persists in the superconducting 
state in the underdoped regime. A theoretical analysis using a new sum-rule reveals that the 
superconducting state cannot be described within the standard BCS-Migdal Eliashberg the- 
ory applied with a pure (cos kx—cos ky) gap. We have also established a new relation between 
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the behavior of the low-energy Raman response and that of the temperature-dependence of 
the superfluid stiffness, both associated with nodal physics. All these observations empha- 
size the key importance of the momentum space differentiation between nodal and antinodal 
regions in both the normal and superconducting state in the proximity to the Mott transi- 
tion. The physics of momentum space differentiation can have multiple origins and has been 
stressed early on by many authors 0, E^- Modern theoretical approaches, such as cluster 
perturbation theory, functional renormalization groups, and cluster extensions of Dynamical 
Mean Field Theory have shown to be able to capture this physics in both the weak coupling 
and strong coupling regimes of the Hubbard model js^, .3^ . Our results provide further 
tests of these methods and should stimulate further development of these approaches, aiming 
in particular at a quantitative evaluation of the ERS response. 



METHODS 



Details of the experimental procedure 

The Hg-1201 single crystals have been grown by the flux method. The detailed procedure 
of the crystal growth is described elsewhere j40|. The magnetically measured transition 
widths are less than 5K for all samples. The ERS experiments have been carried out using 
a triple grating spectrometer (JY-T64000) equipped with a nitrogen cooled CCD detector. 
The crystals were mounted in near-vacuum (10~^mbar) on the cold finger of a liquid helium 
cryostat. Rotation of the crystal was achieved in situ. Red (1.9 eV) and Green (2.4 eV) 
excitation lines have been used for probing the ANR and NR Raman response functions, 
respectively, in order to magnify the ERS signal. It has been experimentally shown that 
resonance effects occur for these lines and are probably due to the matching of the excitation 
energies to the inter-band transitions of Hg-1201 [iH. The spectra presented here have been 
taken below (T = 10 K, corrected from laser heating), and above Tc (T = 100 K for 92K(OD), 
95K(OP) and 89K(UD) crystals, and T = 89 K for 78K(UD) and 63K(UD) crystals). 



Theoretical sum-rule 

The absolute value of the Raman intensity in samples with different doping levels or 
composition is difficult to ascertain experimentally, because it depends on geometric factors 
which are hard to measure, such as the penetration volume of the electromagnetic field into 
the sample. To circumvent this problem, and to make contact with the physics of the Mott 
insulator, we use a sum rule to normalize the Raman intensity. In Ref. ji^], such a sum 
rule was derived for the Hubbard model. This cannot be used to interpret our experiment 
however, because it involves the contribution of the Raman intensity from the upper Hubbard 
band, which cannot be determined experimentally since it occurs at energies where other 
interband transitions take place. For this reason, we have extended the method of Ref. [i^l 
and derived a new low-energy sum rule for the Raman intensity, directly for the t — J 
model (for which the upper Hubbard band has been projected out). This sum-rule relates 
the low energy Raman intensity to the expectation value of a composite operator: 

j dujujx"iuj) = 1<[M, [i/,Mt]] >= E E < X^a'mf,p'U)X,yil)Xss'im) > 

ijkl a-- S' 
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In this expression, Xqc/(2) is the Hubbard operator at site i, and the Greek indices a - ■ - S' 
run over the three possible atomic states 0, t, i- M is the Raman operator (stress tensor) 
M = Ylka kyX„o{'k)Xo^{'k) and H is the t — J Hamiltonain (containing a hopping 

term and a superexchange interaction). We have calculated the tensor T^f^ ^^' appearing 
on the r.h.s of this equation. For a hopping term involving nearest (t) and next-nearest 
neighbor [t'), assuming unbroken spatial and time reversal symmetry, we have checked that 
all non-zero contributions involve at least two fermionic operators (i.e. X^q or Xq^^) living 
on different sites. In the absence of translational symmetry breaking contributions involving 
exactly two such operators are proportional to doping. We have checked this from an 
explicit evaluation using e.g. slave-boson or dynamical mean-field theories. Details will 
be published elsewhere. Hence, the Raman intensity for the t — t' — J model reads: 
o?a;a;x"(<^) = C(t,t' , J) p + ■ ■ ■ at low doping levels, justifying the normalisation used 
above. A similar observation follows from a different sum rule 14311. 
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FIG. 1: Raman responses in the NR (62^) and ANR (Big) as function of doping. The dashed Unes 
on (B2g) (respectively (Bi^)) spectra show the Hnear (resp. cubic) frequency dependences of the 
nodal (resp. antinodal) Raman responses. The arrows indicate the position of the superconducting 
peak maxima. 
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FIG. 2: Universal doping dependence of the ratios hLVAN/ksT^'^^ and Tllon /ksT^'^^ obtained from 
the Big and superconducting Raman peaks. The ratios 2A/fcBrj"°^ deduced from ARPES 
coherent peak in the ANR 
been reported. 
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FIG. 3: Normalized Raman response functions with respect to the sum rule. A weak linear back- 
ground coming from spurious luminescence for intermediate doping, independent of the scattering 
geometry and excitation lines has been subtracted from raw data before performing the normaliza- 
tion (note that without this subtraction the final result is qualitatively similar, i.e the low energy 
slopes of the normalized nodal Raman response are found to be doping independent). 
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FIG. 4: Doping dependence of the slope of the nodal Raman response a = ^(ZA)^, normalized 
to the optimal doping one {p = 0.16). The Fermi liquid par ameter /3 = ^(ZAp)^ extracted form 
the temperature dependence of the penetration depth [liJ, , is also shown, a and /? are both 
found to be doping independent in the range (p = 0.09 — 0.020). 
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FIG. 5: (a) Momentum dependence of the phenomenological superconducting gap function: 
A((/)) = ^ X Ao(i?cos2(/) + {1 — B) cos6(/)). A is related to the maximum amplitude of the super- 
conducting gap, and B to its deviation to the pure d-wave. (b) Momentum dependence of the phe- 
nomenological ZA function: ZA{(f)) = ^/vX x (1 — C cos'^ {2(j))) , where = ^\^=^ = Am(8i? — 6). 
(c) Calculated ANR (Big) response, (d) Calculated NR (62^) response function. Arrows indicate 
the position of the maximum. 
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